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Although the conditions for performing arbitrary unitary operations to simulate the dynamics 
of a closed quantum system are well understood, the same is not true of the more general class of 
quantum operations (also known as superoperators) corresponding to the dynamics of open quantum 
systems. We propose a framework for the generation of Markovian quantum dynamics and study 
the resources needed for universality. For the case of a single qubit, we show that a single nonunitary 
process is necessary and sufflcient to generate all unital Markovian quantum dynamics, whereas a 
set of processes parametrized by one continuous parameter is needed in general. We also obtain 
preliminary results for the unital case in higher dimensions. 



I. INTRODUCTION 

The idea of simulating one quantum system with an- 
other was first suggested in the early 1980's by Manin ||] 
and Feynman A universal quantum computer can 
perform such simulation because of its ability to apply 
arbitrary unitary transformations to arbitrary quantum 
states. The only necessary resources are single qubit 
gates and the controUed-NOT two-qubit gate In fact, 
the controUed-NOT may be replaced by nearly any two- 
qubit interaction Q, and the single-qubit gates can be 
reduced to a finite set Finiteness of the gate set is 
desirable because it reduces the necessary computational 
resources and simplifies the construction of fault tolerant 
gates. 

Using a universal gate set, a quantum computer may 
simulate the time sequence of operations corresponding 
to any unitary dynamics. Such simulation is provably 
efficient and has been implemented in the context of 
NMR quantum computation Q. 

However, quantum systems may undergo interesting 
processes which are not unitary due to interactions with 
their environments. Evolution of such open quantum sys- 
tems is described by quantum operations (or superoper- 
ators). Understanding such dynamics is important for 
studying quantum noise processes js), designing quan- 
tum error correcting codes and performing simula- 
tions of open quantum systems, such as of thermal equi- 
libration 0. 

Clearly, creation of arbitrary quantum operations and 
simulation of arbitrary quantum dynamics using a simple 
set of primitives are desirable goals. However, it is more 
difficult to describe a notion of universality for general 
quantum operations than for unitary operations alone. 
Unlike unitary operations, which form a Lie group, quan- 
tum operations comprise a semigroup due to their irre- 
versibility. The lack of inverse operations for semigroups 



is troublesome, and it is less obvious how best to combine 
quantum operations to form new ones. 

A simple recipe for implementing a general quantum 
operation follows from its unitary representation: any 
quantum operation may be written as a unitary opera- 
tion on an extended system with a trace over the extra 
degrees of freedom. As is well known, this procedure only 
requires an ancillary system of dimension equal to the 
square of the dimension of the system of interest to pro- 
duce arbitrary quantum operations. However, our goal 
is to consider as a resource a small subset of nonunitary 
quantum dynamics applied to the system only, without 
the need to control the extra degrees of freedom. Such re- 
strictions are important in many applications, including 
the experimental simulation of quantum systems. For 
this reason, we exclude the technique arising from the 
unitary representation when building arbitrary nonuni- 
tary quantum dynamics. 

In this paper, we begin to study possible methods for 
simulating the dynamics of open quantum systems by 
some time sequence of operations. We take the approach 
of considering only processes which result from interac- 
tion with a Markovian environment in the Born approxi- 
mation. We refer to this class of dynamics as Markovian 
quantum dynamics and refer to the semigroup they com- 
prise as a Markovian semigroup. Such processes have a 
convenient description in terms of their generators, a con- 
cept analogous to the Hamiltonian of unitary dynamics. 
Therefore, simulation of Markovian quantum dynamics 
is reduced to building generators for Markovian semi- 
groups. 

We define two allowed procedures for transforming 
semigroup generators, linear- combination and unitary 
conjugation. Using these procedures, we show how to 
build more complicated generators from simple ones, and 
we explore in detail the required resources for the case of 
a single quantum bit. 
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The structure of the paper is as follows: in Section || 
we present representations of quantum operations and 
the description of Markovian quantum dynamics in terms 
of semigroups. Then, in Section III, we describe the 



procedures of linear-combination and unitary conjuga- 
tion. Using these procedures, we present universal sets 
of generators for unital Markovian quantum dynamics on 
a qubit in Section [V and for general Markovian quantum 
dynamics on a qubit in Section Finally, we conclude 
with some open questions and directions for further in- 
vestigation. 



II. QUANTUM OPERATIONS AND 
MARKOVIAN SEMIGROUPS 

A quantum state is described by a density matrix p 
which is positive semidefinite and has tvp = 1. The most 
general state change of a quantum system, a quantum 
operation, is a linear map £ which is trace preserving 
and completely positive. £ acts on p to produce a state 
£{p). There are many representations for such a map. 
The operator sum representation 



£{p) = AkpA 



(1) 



(note that we use Einstein summation where appropri- 
ate) — and its corresponding fixed-basis form — is 
convenient in that the constraints of trace preservation 
and complete positivity may be simply expressed. For 
example, complete positivity is inherent in Eq. (|l|) and 
trace preservation is equivalent to = I, where / 

is the identity matrix. However, the composition of two 
operator sum representations is complicated, usually re- 
sulting in a rapidly increasing number of terms. On the 
other hand, a manifestly linear representation 



(2) 



(where M is a matrix with composite indices) makes 
composition of operations trivial, yet obfuscates the con- 
straints @. 

Instead of considering all possible dynamics, we will 
simplify the problem by focusing on Markovian quantum 
dynamics. We describe these processes informally here, 
saving a more complete presentation based on [p^^p^ for 
Appendix Every such process corresponds to some 
interaction which, if applied for a duration t, induces 
a quantum operation £t. The class of quantum opera- 
tions £t forms a Markovian semigroup. The time t may 
vary continuously. The operations must be stationary 
and Markovian, such that 



£s£f — £s+t ■ 



(3) 



Here £s£t denotes composition of the operations, i.e., 
£s o £i. Each Markovian semigroup describes the dy- 
namics resulting from some interaction with a Markovian 
environment in the Born approximation. 



Note that this terminology differs slightly from that 
used elsewhere. For example, Davies does not include the 
constraint of trace preservation when defining a Markov 
semigroup i n |14| and, curiously, uses "Markov" to mean 
"unital" in |15|]. For a precise definition of Markovian 
semigroups as used in this paper, see Appendix ^ 

The advantage of considering only Markovian semi- 
groups is that they are uniquely determined by their gen- 
erators. The generator Z of £t is defined by its action on 
an arbitrary input p, 



Zip) 



lim 

no 



£t{p) 



t 



(4) 



In a sense, Z can be thought of as the "Hamiltonian" 
corresponding to £t ■ Exponentiation gives 



£t = e 



zt 



lim \I--Z 

n 



(5) 



where X is the identity quantum operation. The genera- 
tor also satisfies the differential equation 



dpjt) 
dt 



Z{p{t)). 



(6) 



which is known as a master equation. Through this anal- 
ysis, simulating £t ^t > Q is reduced to simulating its 
generator. 

Gorini, Kossakowski, and Sudarshan have shown that 
Z is the generator of a Markovian semigroup on an N 
dimensional Hilbert space if and only if it can be written 
in the form ||l^ 

Z{p) = ~i[H,p\ + a^^{[F^p,Fl] + [F„,pFj]) , (7) 

where a^p is an {N"^ ^ 1) x (^^ ^ 1) positive matrix (with 
a, /? € [1, iV^ — 1]) and {F^} is a linear basis of traceless 
operators on the density matrices. We refer to the matrix 
aafj as the GKS matrix. For related formulations, such 
as the "diagonal" form introduced by Lindblad (which 
also applies to countably infinite-dimensional systems), 
see fl^Jl^ . Physically, H corresponds to unitary dynam- 
ics which can be produced by a system Hamiltonian as 
well as unitary dynamics induced by a coupling between 
the system and the bath — the so-called Lamb shift. 

It will greatly simplify the discussion to choose a Her- 
mitian basis which is orthonormal under the trace norm. 
Such a basis is assumed for the rest of the paper. There- 
fore, 



tr(F„i^^) 



(8) 



and tr(FQ,) = 0. Note that we can always reduce a GKS 
matrix which is expressed in an overcomplete or non- 
orthonormal traceless basis to a representation involving 
a linearly independent orthonormal traceless basis. 

There are other ways to describe the generator of a 
Markovian semigroup. For example, Z[p) may always 
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be written as an afRne transformation of p, just as any 
quantum operation can be written as a linear transfor- 
mation as in Eq. (||). In this paper, we find it simplest to 
represent generators by the GKS matrix, and we describe 
the relationship between the GKS matrix and the affine 
representation in Appendix |^. 

To make our description of Markovian quantum dy- 
namics concrete, we present some important examples of 
qubit noise processes We choose the basis {F^} to 
be the normalized Pauli operators -^{ca,, cTj,, cTz}, and 
we write the density matrix of a qubit as 



P 



Poo Poi 



\PlO Pll 

The first process, phase damping, acts on a qubit as 



(9) 



\ e ^Vio Pll 



(10) 



where 7 is a decay constant and t is the duration of the 
process. The generator has a GKS matrix with a™ = 
^ and all other = 0. The second example is the 
depolarizing channel, which acts on a qubit as 



l+e "''(poo-pii) 



l-l-e-^'(pii-poo) 



(11) 



Its GKS matrix has the nonzero elements af^^ — 
7/4. Our final example is amplitude 



-,DEP _ „DEP 



'■22 



*33 



damping, which acts on a qubit as 

^AD(^)^^Poo + a^-e-")pn ^-JJ'PoA, (12) 
\ e ' pio e Pll J 

The GKS matrix is given by 



r 



'1 -i 
i 1 
.0 0. 



(13) 



Note that the GKS matrix is real and diagonal for phase 
damping and the depolarizing channel and has rank one 
for phase damping and amplitude damping. 



III. COMPOSITION FRAMEWORK: 
LINEAR-COMBINATION AND UNITARY 
CONJUGATION 

Recall that our goal is to find a simple way of combin- 
ing as few primitive as possible to produce all possible 
£ via some time sequence of operations. To make this 
problem well-posed, we must choose reasonable methods 
for composing quantum operations to make new ones. 
We have not found a simple way to express the compo- 
sition of two semigroup processes of finite duration, and 
such composition need not preserve Markovity. However, 



a natural way to combine semigroup processes is by a pro- 
cedure we call linear- combination: the processes act one 
after another for small amounts of time. In the limit of 
infinitesimal time steps, two processes £^ and £^ can be 
combined to produce 



£r'^ hm (£,%£,"/„)", 



(14) 



where £t~^^ forms a Markovian semigroup if £f and £f 
do. Moreover, if f ° and £f have generators A and B, 
then applying Lie's product formula to the generators. 



lim ( 



(15) 



In other words, the generator of a process formed by 
linear-combination is the sum of the constituent genera- 
tors. The generalization to produce a positive sum of any 
finite number of generators is straightforward. When all 
generators are expressed in the form of Eq. ^ using the 
same basis {Fa} (as we assume for the rest of the paper), 
linear-combination corresponds to a positive sum of the 
GKS matrices of the constituent generators. 

We also assume the capability to apply arbitrary uni- 
tary operations to the system, since these tasks are fea- 
sible and well understood. Using linear-combination, we 
may produce the two terms in Eq. (Q) separately. Assum- 
ing the ability to create any unitary dynamics, it remains 
to generate the second term under the assumption H = 0. 

We now turn to the second procedure to transform the 
GKS matrix, called unitary conjugation. This procedure 
transforms £ according to 



U^£U, 



(16) 



where U{p) — UpU^ for some unitary operator U. Note 
that unitary conjugation preserves all the Markovian 
semigroup properties. We will see that the effect of uni- 
tary conjugation is to apply £ in a different basis, pro- 
ducing a new operation which may be used on its own 
or in linear-combination. To understand how the GKS 
matrix transforms, we prove the following theorem: 

Theorem 1 For an N dimensional system, 
unitary conjugation by U € SU (TV) results 
in conjugation of the GKS matrix by a cor- 
responding element in the adjoint representa- 
tion of SU{N). 

Proof: Suppose the Markovian semigroup has genera- 
tor A and GKS matrix Oa/i. Conjugation by U results in 
the evolution 

U^e^'U = lim (x- ^a\ U (17) 



lim {l--U^Au] 
n I 



(18) 



In other words, the new generator is A' = U'^ ALA. Ex- 
pressed in the form of Eq. (^ (with = 0), we find 
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A'{p)=a^p{[U^F^Up,U^Flu] 
+ [U^F^U,pU^Flu]). 



(19) 



Evidently this unitary conjugation induces a change of 
basis Fa WFaU, which is still Hermitian, orthonor- 
mal, and traceless. We can expand the new basis in terms 
of the old one: 



This implies 



U^FaFlu : 



(20) 



(21) 



Taking the trace of Eq. (|T|) , and using the orthonormal- 
ity ofthef^^ (Eq. (|)), 



(22) 



In other words, Cg-y is a unitary matrix. Further, by 
substituting Eq. ( [2C| ) into Eq. ([l9|), we obtain the trans- 
formed GKS matrix 



(23) 



Denoting the matrices and Cap by A and C, A' — 
AC* . The effect of unitary conjugation is to conjugate 
the original GKS matrix by . 

Note that C is not arbitrary, but is determined by U 
in the following manner. Suppose we choose {Fa} to be 
the generators of SU{N). Then we have 



[Fa,F0] = ifaPfFj , 



(24) 



where fa/3j are the real structure constants for the Lie 
algebra generated by {Fa}- Setting U = e^^'^^'f , we find 
to first order in an infinitesimal that 

U^^FaU = (/ - ir^F^)Fa{I + ir^F^) 



= Fa 



ir^iifjapFfs) 



Thus 



(C'^)q/3 = Sa0 + ir^iif^ap) 



(25) 



(26) 



is in a Lie group generated by the matrices (G'^)^/? = 

It is an elementary fact of group theory that these Gj 
are the generators of the adjoint representation of the 
Fa algebra. Therefore, U = e^^''^'' e SU{N) induces 



conjugation of the GKS matrix by C 



the 



adjoint representation of SU{N). □ 
As an example of an application of these two methods, 
we can perform linear-combination of amplitude damp- 
ing, (Eq. (12)), and damping in the opposite direc- 
tion, X ^ where X{p) = (T^pa^, to simulate gener- 
alized amplitude damping with an arbitrary mixture of 
the ground and excited states as the fixed point. 



IV. UNITAL MARKOVIAN QUANTUM 
DYNAMICS 

We now use the resources we have defined to simulate 
Markovian quantum dynamics. We first consider unital 
processes, those that fix the identity. Well-known uni- 
tal processes on a qubit include phase damping and the 
depolarizing channel. 

We first characterize unitality for GKS matrices. Z 
is the generator of a unital Markovian semigroup iff 
Z{I) = 0. Using Eq. @i (with H = 0), Eq. (||), and 
the antisymmetry of faf3-y, we find 

Z{I) = 2aa0[Fa,Fl] = 2aa0[Fa,Fi3] 

~ "^iaapfap^Fj = 2i ^ (Oq,^ — apa)fal3-tF^ 

Q</3 



-4 ^ lu\(aap)fai3-iF^ 



(27) 



a<l3 



By orthonormality of the Fa, Z{I) =0 iff 

^ Im(a„;3)/„;3^ = V7. (28) 

a<(3 

In general, reality of the GKS matrix is a sufficient con- 
dition for unitality. When N — 2, the sum has only one 
term, so that this condition is also necessary. Thus the 
unital Markovian semigroups on a qubit are exactly those 
generated by real GKS matrices. 

For N > 2, reality of the GKS matrix is not nec- 
essary for the corresponding process to be unital. In 
Appendix we give an example of a unital operation 
of dimension = 3 for which the GKS matrix is not 
real. Since SU{N) for A^ > 3 contains an isomorphic 
copy of SU{3), it follows that there are unital Marko- 
vian semigroups generated by complex GKS matrices for 
all A^ > 3. The set of Markovian semigroups generated 
by real GKS matrices is a proper subset of the set of all 
unital Markovian semigroups. 

We now focus on Markovian quantum dynamics on 
a qubit and consider the effects of unitary conjugation. 
We take -Fi,2,3 = o-x,y,z as before. The generators are 
represented by real positive semidefinite GKS matrices 
A. The transformations induced by unitary conjugation 
are simply 



where 



A' 



G2 



G3 



Ae 



(29) 




(30) 
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can be found using fa[3-y = e-aji-y for the Pauli matrices. 
Note that the G^'s are simply the generators of 5*0(3) 
— as is well known, SO(?)) is the adjoint representation 
of SU{2). 

Having characterized unital Markovian quantum dy- 
namics on a qubit, we are ready to state the following: 

Theorem 2 All unital Markovian quantum 
dynamics on a single qubit can he simulated 
with linear- combination and unitary conjuga- 
tion given phase damping as a primitive op- 
eration. 

Proof: When A is positive semidefinite and real, A = 
ODO^ for some O S 5*0(3) and D diagonal with diago- 
nal entries di > 0. Let A^^ denote the GKS matrix for 
phase damping. Then 

A = ^(diOe^f«MP°e-'^^^^0^ 
7 

+d20e't«MP°e-^^«iO^ + daO^P^O^) . (31) 

Thus the GKS matrix for any unital Markovian semi- 
group can be simulated using linear-combination of the 
unitary conjugates of phase damping. □ 
Note that any unital Markovian semigroup with rank 
one GKS matrix can be used in place of phase damping. 

V. GENERAL MARKOVIAN QUANTUM 
DYNAMICS ON A QUBIT 

We now consider general Markovian quantum dynam- 
ics on a single qubit. We prove the following: 

Theorem 3 To simulate arbitrary Markovian 
quantum dynamics on a qubit with linear- 
combination and conjugation by unitaries, it 
is necessary and sufficient to he able to per- 
form all operations with GKS matrix A{6) = 
d{e)d{ey, where a{e) = (cos 6*, i sin 61,0)^ 
andee [0, f]. 

Sufficiency: Let A be the GKS matrix to be simulated. 
Since A is positive semidefinite, we apply the spectral 
decomposition to express ^4 as a positive sum of outer 
products, 

A = Y, ^fe«fe4 . (32) 
fe 

where Afc > and a's are normalized to unit length. By 
linear-combination, it suffices to simulate all aa^ . Sepa- 
rating into real and imaginary parts, we may write any 
such vector as 

a = -I- ia' . (33) 

Because a only appears in outer products, the overall 
phase of a is irrelevant. In other words, we have a'a'^ — 
da) , where 



a! = e*'^a = {a^ cos — a^ sin 

-\-i{d^ sin Tp -\- d^ cos ip) . (34) 

This transformation maps the two parameters 

/fci = |a«|2-|ci^|2 (35) 
k2 = 2d^^d^ (36) 

according to 

( K \ _ ( cos 21/; -sin2'0\ ( ki\ , . 

\k'2)~\sm2i} cos2i/' J\k2j' ^ ' 

Because we may choose tp arbitrarily, we make the choice 

tan2?/> = -fca/Zci , (38) 

such that fcj = 0; in which case a'^ and a'^ are orthog- 
onal. Moreover, we can choose k[ — ki/ cos2ip > such 
that \d ^\ > \d^\. Thus, without loss of generality, we 
may assume that a has a real part no shorter than its 
imaginary part, and that the two parts are orthogonal. 

Performing a unitary transformation on the operation 
effects conjugation by G G 50(3): 

aa^ Gad^G'^ ^ {Ga){Ga)^ . (39) 

Because G is real, it does not mix the real and imaginary 
parts of a in Ga. In other words, G simultaneously ro- 
tates the two vectors and a' . Therefore, G can always 
be chosen to align with the -\-x axis and with the 
-\-y axis, so that we can write 

Ga = a{e) = (cos 6*, i sin 0, Of , (40) 

where 9 e [0,f]. 

Necessity: The A(0), being rank one, are extreme in 
the convex cone of all the positive matrices. Thus linear- 
combination cannot be used to compose a new A(9'). 

Because scalar multiplication of a by a phase commutes 
with an 50(3) transformation, it suffices to show that 
given a phase and rotation G such that Ge'^^d{9) = 
d{0') with 0,9' e [0,f], then 9 = 9'. To see this, note 
that given A;2 = fcj = and fci , fc^ > in Eq. (|7|) , the 
phase transformation must be trivial. Thus the real and 
imaginary parts are unchanged by e**^, and have to be 
unchanged by G to remain aligned with the -\-x and -\-y 
axes. Therefore, 9' = 9. □ 

Note that the set of required operations includes ampli- 
tude damping, with = a(7r/4) — (1, i, 0)-^/\/2, and 
phase damping about the x axis, with flPD' - a(0) = 
(1,0,0)^. 

We have seen that a one parameter set of generators is 
necessary and sufRcient to simulate all Markovian quan- 
tum dynamics on a qubit. With this result in mind, we 
can estimate the number of parameters in the universal 
set of generators needed for an N dimensional quantum 
system. Generating the (7V^ - 1) x (7V^ - 1) GKS ma- 
trix can be reduced to generating all those with rank one 
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by linear-combination. Thus we have to obtain all (nor- 
malized) complex (7V^ — l)-vectors a up to a phase, with 
2{N'^ — 1) — 2 free parameters. As there are iV^ — 1 degrees 
of freedom in SU{N), unitary conjugation will eliminate 
at most N'^ — 1 parameters, leaving iV^ — 3 parameters to 
be obtained, perhaps by a continuous set of primitives. 
In the case of a qubit (N = 2), we have seen that this 
bound is tight. 

VI. CONCLUSIONS AND OPEN QUESTIONS 

In the search for a way to simulate the dynamics of 
open quantum systems using a simple set of primitives, 
we have set up a framework to study the notion of univer- 
sality for Markovian quantum dynamics. We have shown 
how the generators of Markovian semigroups transform 
under the composition procedures of linear-combination 
(Eq. ([Tsl)) and unitary conjugation (Theorem 1). For the 
case of a single qubit, we have shown that one primitive 
unital operation suffices to simulate all other unital op- 
erations (Theorem 2), and we have exhibited a necessary 
and sufficient single parameter universal set for general 
qubit Markovian quantum dynamics (Theorem 3). 

The most immediate open questions are related to the 
corresponding results for higher dimensions. We have 
seen that for iV > 3, unital operations no longer corre- 
spond to real GKS matrices. An interesting open ques- 
tion is whether there is still a finite generating set for 
unital Markovian quantum dynamics. It would also be 
interesting to see how the results for general processes 
scale to TV > 3 dimensions and how many parameters 
are needed to describe the extreme points of the basis 
set. This will require investigation of the specific adjoint 
representation of SU{N) that acts on {N^ - 1) x {N^ - 1) 
GKS matrices. 

Further open questions arise from other possible com- 
position rules for quantum operations. If we lift the re- 
striction to Markovian quantum dynamics, it is reason- 
able to include composition procedures other than linear- 
combination and unitary conjugation. For instance, we 
might allow direct composition of operations. Or we 
might consider implementing adaptive methods of quan- 
tum control, as suggested by Lloyd and Viola We 
might also combine processes probabilistically (either by 
tossing coins or by performing different operations on dif- 
ferent parts of the sample in a bulk quantum computer) , 
giving a convex combination of quantum operations. 
Such combination does not preserve Markovity in gen- 
eral, but it has other attractive properties; for example, 
any unital operation on a qubit can be written as a con- 
vex sum of unitary operations (although this does not 
hold for higher dimensional systems pl|]). The charac- 
terization of the extremal operations on a qubit presented 
in [ p^ can be used to extend this to nonunital operations. 

Finally, we might consider a formulation of the prob- 
lem which requires only that we be able to come arbitrar- 



ily close to a given quantum operation [£2|. This would 
be more in line with the usual notion of universality for 
unitary operations, and might prove fruitful as a way of 
finding smaller basis sets for general quantum operations. 
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APPENDIX A: FORMAL PROPERTIES OF 
MARKOVIAN SEMIGROUPS 

A strongly continuous one-parameter semigroup on a 
complex Banach space B is defined as a family £t of 
bounded linear operators £t : B B parametrized by 
real i > which satisfy 

(a) So = 1; 

(b) £s£t = £s+t, and 

(c) the map {t,p) £t{p) from [0,cx)) x _B to S is 
jointly continuous. 

The generator Z of a strongly continuous one- 
parameter semigroup is determined by 

Z(p) = lim^^M^. (Al) 

The domain of Z, Dom(Z), is defined to be the space for 
which the above limit exists. Dom(2^) is a dense linear 
subspace of B. If p e Dom(Z), 

^^Z{p{t)) Vt>0, (A2) 

and the semigroup is defined by its generator according 
to 

£t = e^* = Um (l - -z\ , (A3) 

n— ►oo y " / 

the inverse being a bounded operator for sufficiently large 
n. 

A one-parameter semigroup is norm continuous if and 
only if the generator is bounded, in which case 
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£t = lim 



n-^oo -* — ' n! 



(A4) 



Continuous one parameter semigroups capture the 
Markovian and stationarity features of the Markovian 
quantum dynamics of interest. The remaining features to 
be incorporated are complete positivity and trace preser- 
vation. This leads to our definition: 

A Markovian semigroup is a norm continuous 
one-parameter semigroup of completely posi- 
tive, trace preserving linear maps. 



Conversely, given {L,p}, Z is uniquely determined. For 
a fixed basis {Fa}, a unique decomposition in the form 
of Eq. (Q) exists, so that the GKS matrix is uniquely de- 
t erm ined. In other words, the linear system Eqs. (B5)- 
(B4) can be inverted. Thus the {L,p} that represent 



Markovian semigroups are those for which a solution 
to Eqs. (p35D-(B4) exists and corresponds to a positive 
semidefinite GKS matrix. 



APPENDIX C: A UNITAL PROCESS WITH 
COMPLEX GKS MATRIX 



APPENDIX B: THE AFFINE REPRESENTATION 
AND THE GKS MATRIX 

Following the discussion in Section ||, consider a basis 
for traceless operators {Fa} that is Hermitian and trace 
orthonormal. Then we can express a density matrix as 
p = po I -\- PaFa where po, pa are real numbers. Due 
to trace preservation, the linear representation of the gen- 
erator 2^ of a Markovian semigroup can be reduced to an 
affine map on the traceless components only: 



Pa = LapPl3 + Pa 



(Bl) 



1) X {N'^ - 1) 
l)-dimensional 



In Eq. (Bl), Laf3 are entries of an {N"^ 
matrix L, and pa are entries of an {N"^ 
vector. 

In the qubit case, let {Fa} = ■^{ax,CFy,az}- Using 
Eq. (Q), we obtain a one-to-one correspondence between 
the GKS matrix A (with entries Oafj) and {L,p} in the 
affine representation: 



-2 (022 + 033) 

ai2 + 021 - 

ai3 + asi 
A + A^ ~{2 tiA) I 



012 + 021 
2 (0114-033) 

^23 + ^32 



Ol3 + O31 
O23 + 032 
-2 (011-I-O22J 



p = 4(Im(o32),Im(ai3),Im(o2i)) . 

As positivity of the GKS matrix is equivalent to complete 
positivity of Eq. (^ , the correspondence between A and 
{L,p} allows complete positivity of the afhne represen- 
tation to be easily characterized. 

For higher dimensions, the affine representation {L,p} 
and the GKS matrix for a generator Z are still in one-to- 
one correspondence. In particular, let [Fa,Fp\ = ifafS-yFj 
and {FajFp} = hap-^F^ where fap-y are the real structure 
constants and haPf are real. Then the entries of L and 
p are given by 



-hajxfxpr] Ini(Oa/3) 
-/q7a/a/3i) Re(0Q/3) , 



A convenient set of generators for SU (3), known as the 
Gell-Mann matrices, are defined as 




V2X2 = 



V2A4 



■Q 0^ 

i 

,0 0. 

1\ 



.1 0/ 

0\ 

1 



0/ 
^ 


,0 -2, 



.0 1 
1 
1 



(CI) 



The only nonvanishing structure constants are /321 = 

— \/2, /651 = I/V2, /642 = /74I = /752 = —^/V^, /854 = 

/876 — ~-V3/V2 together with cyclic permutations of the 
indices. 

The following process defined over the Gell-Mann ma- 
trices is unital, although A is complex: 



A^- 
2 



/O 













Vo 







1 




— z 




0\ 





i 

1 
1 
1 

000/ 



(C2) 



The eigenvalues are 2 and with degeneracies 2 and 6, so 
the process is well defined. To show that it is unital, we 
need to show J2a<f3 '^a/sfap-y = V7. First, Oap = for 
all a < P except for 047 = 050 = i. The criterion reduces 
to /47^ + /se^ = V7. But /47^ = /se^, = V7 7^ 1, 
and /471 — 1/2, /561 = —1/2. Therefore, the process is 
unital. 



(B4) 



Q</3 



/Q/3^Im(Oa^) . 



(B5) 
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